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ABSTRACT: A boundary layer analysis is presented for the effect of radiation on mixed convection flow of a non- 
Newtonian nanofluid in a nonlinearlv si retching sheet with heat source/sink. The inicropohir model is chosen for the non- 
Newtonian fluid since lite spinning motion of the nanoparticles as they move along the streamwise direction can he best 
described by the micropolar fluid model. Numerical results for friction factor, surface heat transfer rate and mass transfer 
rate have been presented for parametric variations of the micropolar material parameters, Brownian motion parameter Ng, 
thermophoresis parameter Nt, radiation parameter R, heat source/sink parameter O and Schmidt number Sc. The 
dependency of the friction factor, surface heat transfer rale and mass transfer rate on these parameter:, has been discussed. 
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I. INTRODUCTION 

Study of non-Newtonian fluids over a stretching surface achieved great attention due to its large number of 
application. In fact, the effects of non-Newtonian behavior can be determined due to its elasticity, but sometimes rheological 
properties of fluid are identified by their constitutive equations. In view of rheological parameters, the constitutive equations 
in the non-Newtonian fluids are morecomplex and thus give rise the equations which are complicated than the Navier-Stokes 
equations. Many of the fluids used in the oil industry and simulate reservoirs are significantly non-Newtonian. In different 
degree, they display shear-dependent of viscosity, thixotropy and elasticity (Pearson and Tardy [1]; Ellahi and Afza [2]; 
Ellahi [3]). Gorla and Kumari [4] studied the mixed convection flow of a non-newtonian nanofluid over a non-linearly 
stretching sheet. 

Nanoparticles range in diameter between 1 and 100 nm. Nanofluids commonly contain up to a 5% volume fraction 
of nanoparticlesto ensure effective heat transfer enhancements. One of the main objectives of using nanofluids is to achieve 
the best thermal properties with the least possible (<1%) volume fraction of nanoparticles in the base fluid. The term 
nanofluid was first proposed by Choi [5] to indicate engineered colloids composed of nanoparticles dispersed in a base fluid. 
The characteristic feature of nanofluids is thermal conductivity enhancement; a phenomenon observed by Masuda et al. 
[6]. There are many studies on the mechanism behind the enhanced heat transfer characteristics using nanofluids. Eldabe et 
al. [7] analyzed the effects of magnetic field and heat generation on viscous flow and heat transfer over a nonlinearly 
stretching surface in a nanofluid. 

Micropolar fluids are subset of the micromorphic fluid theory introduced in a pioneering paper by Eringen[8]. 
Micropolar fluids are those fluids consisting of randomly oriented particles suspended in a viscous medium, which can 
undergo a rotation that can affect the hydrodynamics of the flow, making it a distinctly non-Newtonian fluid. They constitute 
an important branch of non-Newtonian fluid dynamics where microrotation effects as well as microinertia are exhibited. 
Eringen's theory has provided a good model for studying a number of complicated fluids, such as colloidal fluids, polymeric 
fluids and blood. 

In the context of space technology and in processes involving high temperatures, the effects of radiation are of vital 
importance. Studies of free convection flow along a vertical cylinder or horizontal cylinder are important in the field of 
geothermal power generation and drilling operations where the free stream and buoyancy induced fluid velocities are of 
roughly the same order of magnitude. Many researchers such as Arpaci [9], Cess [10], Cheng and Ozisik [11], Raptis [12], 
Hossain and Takhar [13, 14] have investigated the interaction of thermal radiation and free convection for different 
geometries, by considering the flow to be steady. Oahimire et al.[15] studied the analytical solution to mhd micropolar fluid 
flow past a vertical plate in a slip-flow regime in the presence of thermal diffusion and thermal radiation. El-Arabawy [16] 
studied the effect of suction/injection on a micropolar fluid past a continuously moving plate in the presence of radiation. 
Ogulu [17] studied the oscillating plate -temperature flow of a polar fluid past a vertical porous plate in the presence of 
couple stresses and radiation. Rahman and Sattar [18] studied transient convective heat transfer flow of a micropolar fluid 
past a continuously moving vertical porous plate with time dependent suction in the presence of radiation. Mat et al. [19] 
studied the radiation effect on marangoni convection boundary layer flow of a nanofluid. 

The heat source/sink effects in thermal convection, are significant where there may exist a high temperature 
differences between the surface (e.g. space craft body) and the ambient fluid. Heat Generation is also important in the 
context of exothermic or endothermic chemical reaction. Sparrow and Cess [20] provided one of the earliest studies using a 
similarity approach for stagnation point flow with heat source/sink which vary in time. Pop and Soundalgekar [21] studied 
unsteady free convection flow past an infinite plate with constant suction and heat source. Rahman and Sattar [22] studied 
magnetohydrodynamic convective flow of a micropolar fluid past a vertical porous plate in the presence of heat 
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generation/absorption. Lin et al. [23] studied the marangoni convection flow and heat transfer in pseudoplastic non- 
newtonian nanofluids with radiation effects and heat generation or absorption effects. 

The present work has been undertaken in order to analyze the two mixed convection boundary layer flow of a non- 
Newtonian nanofluid on a non-linearly stretching sheet in the presence of radiation and heat source/sink. The micropolar 
model is chosen for the non- Newtonian fluid since the spinning motion of the nanoparticles as they move along the 
streamwise direction can be best described by the micropolar fluid model. The effects of Brownian motion and 
thermophoresis are included for the nanofluid. The governing boundary layer equations have been transformed to a two- 
point boundary value problem in similarity variables and the resultant problem is solved numerically using the fourth order 
Runga-Kutta method along with shooting technique. The effects of various governing parameters on the fluid velocity, 
temperature, concentration, and Nusselt number are shown in figures and analyzed in detail. 

II. MATHEMATICL ANALYSIS 

Consider a steady, mixed convection boundary layer flow of a micropolar nanofluid over a non-linearly stretching 
sheet. The velocity components in x and y directions are u and v respectively. At the surface, the temperature T and the nano 
particle fraction take values T w and C w , respectively. The ambient values, attained as the radial distance r tends to infinity, of 
T and C are denoted by T m and , respectively. The governing equations within boundary layer approximation may be 
written as: 
Continuity equation 
du dv 



8x dy 
Momentum equation 



O (2.1) 



du du u+K)d A U KdN 



I P Jdy 2 P % 
Angular Momentum equation 

Energy equation 



d l r_a rrL dq L+ \ ecerf pr_Y ffr) 



"W^[" B ^\^K^) yipcff^-^ {2A) 

Species equation 

"f + vf=Z^ + f^ (2 . 5) 

The boundary conditions may be written as: 

u=u w =cx n ,v=0,N=-m— ,T=T W ,C=C W at y=0 (2.6) 

dy 

«^0,iV-K),r^r QO ,C^C Q o as .y— xx> (2.7) 
In the above equations, N is the component of microrotation vector normal to the xy-plane; U x the free stream velocity, D B 
is the Brownian diffusion coefficient, D T is the Thermophoretic diffusion coefficient, a m is the thermal diffusivity, q r is 
the heat flux, Q 0 is the heat generation or absorption coefficient such that Q Q > 0 corresponds to heat generation while 
<2o < 0 corresponds to heat absorption, g is the acceleration due to gravity, p f is the material density of the base fluid, p p 
is the material density of nanoparticles, (pc)^ ar, d (pc) are the heat capacity of the base fluid and the effective heat 

capacity of the nano particle material, respectively and jU, K, p, j, y and Ct m are respectively the dynamic viscosity, vortex 
viscosity (or the microrotation viscosity), fluid density, microinertia density, spin gradient viscosity and thermal diffusivity. 
We follow the work of many authors by assuming that y=(ju+K/2)j=ju(l+K/2)j where K=k/ju is the material parameter. 
This assumption is invoked to allow the field of equations to predict the correct behavior in the limiting case when the 
microstructure effects become negligible and the total spin reduces to the angular velocity (see Ahmadil4 or Yiicell5). 

The radiative heat flux q r is described by Roseland approximation such that 

4cr* dT 4 

qr= ~^rr^^> — 

(2.8) 

where cr* and & 'are the Stefan-Boltzmann constant and the mean absorption coefficient, respectively. Following Chamkha 
[20], we assume that the temperature differences within the flow are sufficiently small so that they T can be expressed as a 
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linear function after using Taylor series to expand T about the free stream temperature T x and neglecting higher-order 
terms. This result is the following approximation: 



T*&4T£r-3T£ 
In view of equations (2.8) and (2.9), equation (2.4) reduces to 

t3 i ^2 T 



(2.9) 



dT 5T = (,,160-*^ 
" dx + V dy ~ "n 



Proceeding with the analysis, we de.ne the following transformations: 

f^2u 1 n. ( . h(n+l) (n-l)/2 

n ~ 1)/2 "™> 0 ™= T~-Z> = c'-ZL 

The mass conservation equation (2.1) is satisfied by the Cauchy-Riemann Equations 



M 2u 



dy dx 

The Darcian velocity components: 

(pc)f ' (/3c)y ' * u 2 

Re^^^,^^^^^^--^^,^,^^^^^^--^ 3 (2.13) 
" f (pc)fa m Re 2 (.pc)fa m T ao 

R=^,Q= ®2 ^-,Sc^-,S=^^ 

In view of the Equation (2.11), (2.12) and (2.13) , and the equations (2.2), (2.3), (2.5) and (2.10) reduce to the following non- 
dimensional form 



{ i + K/2)g" + fg<-^rg-—(2 g+ n=o 

n+l n+l 

Sc ScN B 
Also the boundary conditions in (2.5) reduces to 

/(0) = 0, / '(0) = 1, g(0) = -mf "(0), 9(0) = m = 1 

/ '(oo) -> 0, g(oo) -> 0, fl(<») -> 0, *go) -> 0 
The wall shear stress is given by 

T w =(ju+k)( 'g0 +kN y=0 
Friction factor is given by 

Cr = T ? w = 2Re; 1/2 f "i 1 ! [(l+tf)/"(0)+Kg(0)] 
The wall couple stress is given by 

The dimensionless wall couple stress may be written as 
V(0) 



(2.14) 

(2.15) 
(2.16) 
(2.17) 

(2.18) 



Heat transfer rate is given by 
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(2.23) 



Local Sherwood number is given by 




A/2 



010) 



(2.24) 



III. 



SOLUTION OF THE PROBLEM 



The set of coupled non-linear governing boundary layer Equations (2.14) - (2.17) together with the boundary 
conditions (2.18) are solved numerically by using Runge-Kutta fourth order technique along with shooting method. First of 
all, higher order non-linear differential Equations (2.14) - (2.17) are converted into simultaneous linear differential equations 
of first order and they are further transformed into initial value problem by applying the shooting technique (Jain et al. [25]). 
The resultant initial value problem is solved by employing Runge-Kutta fourth order technique. The step size Ar/ = 0.05 is 
used to obtain the numerical solution with five decimal place accuracy as the criterion of convergence. From the process of 
numerical computation, the skin-friction coefficient, the wall couple stress, the Nusselt number and the Sherwood number, 
which are respectively proportional to /"(0), g '(0) , — $'(0) and — 0'(O) , are also sorted out and their numerical values 
are presented in a tabular form. 



The governing equations (2.14) - (2.17) subject to the boundary conditions (2.18) are integrated as described in 
section 3. In order to get a clear insight of the physical problem, the velocity, temperature and concentration have been 
discussed by assigning numerical values to the governing parameters encountered in the problem. 

In order to assess the accuracy of the results, we compared our results for R=Q=N T =N B =0 with literature data from 
Gorla and Kumari [4] in Table I and found that they are in excellent agreement. This suggests that the present results are 
accurate. Tables' II-X present data for friction factor, heat and mass transfer rates for parametric values of variables 
governing the problem. Thermophoresis parameter, N T appears in the thermal and concentration boundary layer equations. 
As we note, it is coupled with temperature function and plays a strong role in determining the diffusion of heat and 
nanoparticle concentration in the boundary layer. Table II shows that as the thermophoresis parameter increases, friction 
factor increase, where as the heat and mass transfer rates decrease. N B is the Brownian motion parameter. Brownian motion 
decelerates the flow in the nanofluid boundary layer. Brownian diffusion promotes heat conduction. The nanoparticles 
increase the surface area for heat transfer. Nanofluid is a two phase fluid where the nanoparticles move randomly and 
increase the energy exchange rates. Brownian motion reduces nanoparticle diffusion. Table III shows that as the Brownian 
motion parameter N B increases, mass transfer rate increase, whereas the friction factor and heat transfer rate decrease. Table 
IV shows that as the power law exponent n increases, heat transfer rate increases whereas friction factor and mass transfer 
rate decrease. Table V shows that as the parameter S increases, friction factor, heat and mass transfer rates increases. Table 
VI shows that as the vortex viscosity parameter K increases, the friction factor, and the heat and mass transfer rates increases. 
Table VII shows that friction factor, heat and mass transfer rates increases with the buoyancy parameter X. Table VIII shows 
that as the Schmidt number increases, the friction factor and mass transfer rates increase. Table IX shows that as the constant 
m increases, the wall friction, heat and mass transfer rates decreases. The value m=0 represents concentrated particle flows in 
which the particle density is sufficiently great that microelements close to the wall are unable to rotate. This condition is also 
called as the strong interaction. When m = 0.5 the particle rotation is equal to fluid viscosity at the boundary for one particle 
suspension. When m = 1, we have flows which are representative to turbulent boundary layers. Table X shows that as the 
radiation and heat source/sink parameters increases, the friction factor and mass transfer rate increase, where as heat transfer 
rate decreases. 

Figures l(a)-l(d) and 2(a)-2(d) show the effect of the thermophoresis number, N T as well as the Brownian motion 
parameter, N B on velocity, angular velocity, temperature and concentration profiles. As N T increases, velocity, angular 
velocity, temperature and concentration within the boundary layer increase. As N B increases, temperature increases whereas 
the velocity, angular velocity and concentration decrease within the boundary layer. Figure 3(a)-3(d) shows that as the 
exponent n increases, velocity decreases whereas the angular velocity, temperature and concentration increase. The effects of 
the parameter S on the boundary layer profiles are opposite to the effects of n as seen from Figure 4(a)-4(d). Figures 5(a)- 
5(d) and 6(a)-6(d) show that as the vortex viscosity parameter K increases, the velocity and angular velocity increases 
whereas the temperature and concentration within the boundary layer decrease. As the mixed convection parameter X 
increases, the velocity increases whereas the angular velocity, temperature and concentration within the boundary layer 
decrease. Figure 7(a)-7(d) shows that as the Prandtl number Pr increases, velocity and temperature decreases whereas the 
angular velocity and concentration increase. As the Schmidt number increases, the velocity and concentration boundary layer 
thickness decreases whereas angular velocity and temperature increase. This is seen from Figure 8(a)-8(d). As the constant m 
increases, the velocity decreases whereas the angular velocity, temperatures and concentration increases. This is seen from 
Figure 9(a)-9(d). Figures 10(a)-10(d) and 1 l(a)-l 1(d) shows that as the radiation parameter R and heat source/sink parameter 
Q increases, angular velocity decreases whereas velocity, temperature and concentration increase. 



IV. 



RESULTS AND DISCUSSIONS 
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Table I: Numerical values of (1 + K 1 2)f "(0) and — 0 '(0) at the sheet for different values of K and n when 
A, = 1 .0, S = N T = N B =R = Q = 0.0, m = 0.5 and Pr=0. 72, Comparison of the present results with that of Gorla and 



Kumari [4] 



K 


n 


Present results 


Gorla and Kumari [4] 


(l + K/2)f"(0) 


-0'(O) 


(l+K/2)f"(0) 


-0'(O) 


0.0 


0.50 


-0.129346 


0.586095 


-0.129016 


0.584934 


0.0 


0.75 


-0.312889 


0.568677 


-0.311993 


0.567454 


0.0 


1.00 


-0.444373 


0.555571 


-0.443205 


0.554359 


0.0 


1.50 


-0.621160 


0.536891 


-0.619938 


0.535696 


0.0 


3.00 


-0.874189 


0.507159 


-0.873003 


0.505097 


1.0 


0.50 


-0.293989 


0.593002 


-0.292302 


0.592127 


1.0 


0.75 


-0.497473 


0.577915 


-0.495675 


0.577132 


1.0 


1.00 


-0.643706 


0.566632 


-0.641769 


0.565834 


1.0 


1.50 


-0.840768 


0.550693 


-0.838697 


0.549797 


1.0 


3.00 


-1.123040 


0.525838 


-1.121011 


0.524109 



Table II: Numerical values of / "(0), g '(0), —0 '(0) and — (f> '(0) at the sheet for different values of N T when 

K = S = A = 1,N B =0.3,7M = 0.5,2 = 0.2,« = 2,Pr = 10 andSc=l0. 



N T 


/"(0) 


g\0) 


-0X0) 


-f(0) 


0.1 


-0.664881 


-0.243793 


1.065830 


2.26034 


0.2 


-0.655741 


-0.239314 


1.045970 


2.16199 


0.3 


-0.646809 


-0.234939 


1.026430 


2.07338 


0.4 


-0.638081 


-0.230668 


1.007200 


1.99418 


0.5 


-0.629554 


-0.226498 


0.988283 


1.92407 



Table III: Numerical values of / "(0), g '(0),— 0'(O) and—(/)\0) at the sheet for different values of N B when 



K = S = A = l,N T = 0.3,m = 0.5,(? = 0.2,rc = 2,Pr = 10 andSc=10. 



N B 


/"(0) 


g\0) 


-9X0) 


-f(0) 


0.1 


-0.599733 


-0.211426 


1.111610 


1.30727 


0.2 


-0.635423 


-0.229303 


1.065070 


1.88328 


0.3 


-0.646809 


-0.234939 


1.026430 


2.07338 


0.4 


-0.652032 


-0.237504 


0.990473 


2.16771 


0.5 


-0.654779 


-0.238841 


0.956098 


2.22387 



Table IV: Numerical values of / "(0), g '(0), —0X0) and — (f> '(0) at the sheet for different values of n when 

K = S = R = A = 1,N T = N B = 0.3, m = 0.5, Q = 0.2, Pr = 10 and Sc=10. 



n 


f'XO) 


gXO) 


-0X0) 


-0X0) 


0.5 


-0.197829 


0.172956 


1.00874 


2.20583 


1.0 


-0.430952 


-0.0263596 


1.06582 


2.10284 


1.5 


-0.563788 


-0.15134 


1.10126 


2.04072 


2.0 


-0.649813 


-0.23634 


1.12576 


1.99916 


3.0 


-0.754761 


-0.344031 


1.15695 


1.94704 


5.0 


-0.857058 


-0.452846 


1.18866 


1.89471 



Table V: Numerical values of / "(0), g '(0), — 0 '(0) and — (f) '(0) at the sheet for different values of S when 
K = A = R = l,N T =N B =0.3,m = 0.5,(2 = 0.2,n = 2,Pr = 10 andSc=10. 



s 


f'XO) 


gXO) 


-0X0) 


-f(0) 


0.5 


-0.705108 


-0.267333 


1.11713 


1.99127 


0.6 


-0.694000 


-0.261128 


1.11888 


1.99287 


0.7 


-0.682917 


-0.254928 


1.12062 


1.99445 


0.8 


-0.671858 


-0.248728 


1.12234 


1.99603 


0.9 


-0.660824 


-0.242532 


1.12406 


1.99760 


1.0 


-0.649813 


-0.230340 


1.12576 


1.99916 
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Table VI: Numerical values of / "(0), g '(0), —0'(O) and — (/) '(0) at the sheet for different values of K when 

S = A = R = l,N T =N B = 0.3,m = 0.5,<2 = 0.2,n = 2,Pr = 10 andSc=lO. 



K 


/"(0) 


*'(0) 


-6X0) 


-w 


0.0 


-0.708256 


-0.368713 


0.99671 


2.05477 


0.5 


-0.678646 


-0.284347 


1.01363 


2.06507 


1.0 


-0.646809 


-0.234939 


1.02043 


2.07338 


2.0 


-0.593866 


-0.176569 


1.04481 


2.08581 


3.0 


-0.554020 


-0.142405 


1.05753 


2.09470 



Table VII: Numerical values of / "(0), g '(0), — 0 '(0) and — (f) '(0) at the sheet for different values of X when 

K = S = R = l,N T =Ar B =0.3,m = 0.5,<2 = 0.2,« = 2,Pr = 10 andSc=l0. 



X 


/"(0) 


g'(0) 


-0X0) 


-0X0) 


0.0 


-0.902584 


-0.3743770 


0.976165 


2.03765 


0.5 


-0.772424 


-0.3042750 


1.002770 


2.05617 


1.0 


-0.646809 


-0.2349390 


1.026430 


2.07338 


1.5 


-0.524913 


-0.1662230 


1.047850 


2.08950 


2.0 


-0.406157 


-0.0980246 


1.087500 


2.10473 



Table VIII: Numerical values of / "(0), g '(0), — 0 '(0) and — (f) '(0) at the sheet for different values of X when 
K = S = R = l,N T =TV B =0.3,m = 0.5,G = 0.2,n = 2,Pr = 10 andSc=10. 



Sc 


f'XO) 


gXO) 


-0X0) 


-AO) 


1 


-0.444950 


-0.144039 


1.27846 


-0.185624 


10 


-0.649813 


-0.236340 


1.12596 


1.999160 


20 


-0.683747 


-0.254022 


1.10025 


3.160030 


50 


-0.714107 


-0.270606 


1.08001 


5.353600 


100 


-0.729246 


-0.279186 


1.07134 


7.757390 



Table IX: Numerical values of / "(0), g '(0), — 0 '(0) and — (f) '(0) at the sheet for different values of m when 



K = S = A = R = l,N T =N B = 0.3,(2 = 0.2,n = 2,Pr = 10 and Sc=10. 



m 


f'XO) 


gXO) 


-0X0) 


-f(0) 


0.0 


-0.556371 


0.175433 


1.14406 


2.01556 


0.5 


-0.649813 


-0.23634 


1.12576 


1.99916 


1.0 


-0.781214 


-0.805276 


1.09812 


1.97507 



Table X: Numerical values of / "(0), g '(0), — 0 '(0) and — (f> '(0) at the sheet for different values of R and Q when 



K = S = A = 1,N T =N B =0.3,m = 0.5,« = 2,Pr = 10 andSc=10. 



R 


Q 


f'XO) 


gXO) 


-0X0) 


-^'(0) 


0 


0.2 


-0.677021 


-0.248789 


1.420510 


1.885870 


1 


0.2 


-0.646809 


-0.234939 


1.026430 


2.073380 


2 


0.2 


-0.624100 


-0.225284 


0.825949 


2.165100 


3 


0.2 


-0.605761 


-0.217902 


0.701209 


2.219370 


1 


-0.5 


-0.664790 


-0.243369 


1.701410 


1.556840 


1 


0.0 


-0.655018 


-0.238772 


1.308890 


1.308890 


1 


0.5 


-0.639733 


-0.231650 


0.807817 


2.234470 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

Fig.l(a) Velocity profiles for different values of N T 



9 




n 

Fig.l (b) Angular velocity profiles for different values of N T 




n 

Fig. 1(c) Temperature profiles for different values of N T 























Tl 




Fig.l(d) Concentration profiles for different values of N T 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 



Fig.2(a) Velocity profiles for different values of N B 
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Fig.2 (b) Angular velocity profiles for different values of N B 




Fig.2(c) Temperature profiles for different values of N B 
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Fig.3(c) Temperature profiles for different values of n 
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Fig.3 (d) Concentration profiles for different values of n 
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Fig.4 (b) Angular velocity profiles for different values of S 
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Fig.5 (c) Temperature profiles for different values of K 
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Fig.5(d) Concentration profiles for different values of K 
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Fig.6 (a) Velocity profiles for different values of X 
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Fig.6(b) Angular velocity profiles for different values of X 
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Fig.6(c) Temperature profiles for different values of X 
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Fig.6(d) Concentration profiles for different values of X 
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Fig.7(d) Concentration profiles for different values of Pr 
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Fig.8(a) Velocity profiles for different values of Sc 
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Fig.8(b) Angular velocity profiles for different values of Sc 
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Fig.8(c) Temperature profiles for different values of Sc 
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Fig.8(d) Concentration profiles for different values of Sc 
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Fig.9(a) Velocity profiles for different values of m 
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Fig. 11(c) Temperature profiles for different values of Q 
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Fig. 11 (d) Concentration profiles for different values of Q 

V. CONCLUSIONS 

In this paper, we have presented a boundary layer analysis for the mixed convection flow of a non-Newtonian 
nanofluid on a non-linearly stretching sheet by taking radiation and heat source/sink into account. The micropolar model is 
chosen for the non-Newtonian fluid since the spinning motion of the nanoparticles as they move along the streamwise 
direction can be best described by the micropolar fluid model. The governing boundary layer equations have been 
transformed to a two-point boundary value problem in similarity variables and the resultant problem is solved numerically 
using the fourth order Runga-Kutta method along with shooting technique. The particular solutions reported in this paper 
were validated by comparing with solutions existing in the previously published paper. Our results show a good agreement 
with the existing work in the literature. The results are summarized as follows 

1. The thermophoresis parameter increases, friction factor increase, where as the heat and mass transfer rates decreases. 

2. Brownian motion reduces nanoparticle diffusion. 

3. The radiation and heat source/sink increases the friction factor and mass transfer rate and reduces the heat transfer rate. 

4. The radiation and heat source/sink reduces the angular velocity, and increase the velocity, temperature and 
concentration. 
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